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APPENDIX I 



% Creating a network topology 

network Jopo = topoCinit'); 
addlink(network Jopo) ; 
labelnames(network_topo) ; 

save network__topo; 



% graphically place nodes on screen 
% graphically connect up nodes 
% graphically label nodes 

% save network_topo for future use 



% Top level procedure to compute paths that optimize use of network capacity 
% inputs: 

% D = traffic demand matrix 

% (retrieved from predictions stored in TMS Statistics Repository) 
% network Jopo topology object defining the network topology 
% P = network policy information 

% (matrix of reserved capacity, which indicates links whose use 

% is administratively prohibited or which should not be 

% completely allocated) 

% outputs: 

% allocated_j)aths() = list of paths to set up, to TMS signalling system 



C = capacity (networkjopo); % retrieve network topology information 

C = C - P; 

saved^C = []; 
saved_SLA== []; 
assigned_paths = []; 
round = 0; 

[SLA, S] = create_ordered_sla(D); 

F-SLA(l) 

for F = SLA', 

round = round +1; 
saved_C{ round} = C; 
saved_SLA{round} ^ F; 

F % display the flow 

W = calc_weights('calcweight2^F5C); 
[dist, P] = floyd(W); 
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path = fmdpath(P,F,i,F.j); 

assigned_paths {round}, path = path; 
assigned_paths {round}. flow = F; 

if (isempty(path)) 

fprintf(i;no path for flow:\n'); F 

else 

C = compute_residual_capacity('c - F.bw\path,F,C); 

end 

end 

function [W] = calc_weights(func,F,C) 
% function [W] = calc_weights(func,F,C) 

% 

% Compute the weights by calUng func on each elt of C 

% func must be of the form double func(Flow F, Capacity_eh c, node i, node j) 

func = fcnchk(func); 

fori = l:size(C,l) 

for j l:size(C,2) 

W(ij) = feval(func,F,C(ij),id); 

end 

end 

function [w] = calcweight2(F,c,ij) 
% function [w] = calcweight2(F,c4j) 
% basic weight calc 

if(0 = c) 

w = inf; 
return; 

end 



% rule out paths that can't hack it 
if (F.bw>c) 

w = inf; 

return; 

end 
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^ = 1 / - F.bw) ; % fill links with most capacity first 

function [C] = compute_residual_capacity(func, path, F, C) 
% function [C] = compute__residual_capacity(func, path, F, C) 
% 

% Update capacity characteristics in C to reflect flow F being 

% allocated along path using function func 

% func should be of the form 

% C__element func(C_element c, Flow F) 



if(length(path) 1) 
return; 

end 

func ^ fcnchk(func,'cVF'); 

index = 1 ; 

src = path(index); 

index = index + 1 ; 

for index = index :length(path) 
dst = path(index); 

C(src,dst) - feval(func,C(src,dst),F); 
src = dst; 

end 

function [SLA, S] = create_ordered_sla(D) 

% function [SLA] = create_ordered_sla(D) 

% takes the demand matrix and returns a list of SLAs, 

% SLA of the form [ struct ; struct ; ... ] where struct is [BW, i, j] 

% S of the form [ [BW, i, j] ; [BW, i, j] ; ...] 



s = []; 

for i = l:size(D,l) 

for j = l:size(D,2) 

if(D(i,j) 0) 

S = [[D(iJ)ij];S]; 

end 

end 

end 
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[Y, I] == sortrows(S,l); 

S = Y(size( Y, 1 ) : - 1 : 1 , :); % reverse order 

5 

SLA = struct(W,num2cell(S(:,l));i',nurn2cell(S(:,2));j',num2cell(S(:,3))); 
return; 

10 flmction [path] = findpath(P,i,j) 

% function [path] = findpath(P) 

% 
% 

15 path=[]; 

ifa==j) 

path= [i]; 
return; 

20 end 

if(0 = P(iJ)) 
path=[]; 

else 

25 path=[findpath(P,i,P(i,j))j]; 
end 

function [D, P] = floyd(W) 
% function [D, P] = floyd(W) 
30 % given weights Wij, compute min dist Dij between node i to j 

% on shortest path from i to j, j has immeadiate predecessor Pij 

n = size(W,l); 
if(n~=size(W,2)) 
3 5 error('Input W is not square?? ! ! ') ; 

end 

D = W; 

40 P = repmat([l:n]',[l n]); 

P = P .* ~isinf(W); 
P = P .* ~eye(n); 

for k = 1 :n 
45 for i = 1 :n 
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forj = l:n 

altjjath = D(i,k) + D(k,j); 
if(D(io)>alt_path) 

D(i,j) = alt_path; 

P(ij) = P(k,j); 

end 



end 

end 

k; 

10 D; 

P; 

end 
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APPENDIX II 

function addlink(TOPO) 
% addlink(TOPO) 
% 

% interactively add links to the TOPO 



update(TOPO); 
c_src = 1 ; 
c_dst = 2; 
10 cbw = 3; 



figure{TOPO .cur_fig) 
while- (1) 

fprintf(l,'\n\nHit Button 3 to end...\n\n'); 



% find coords and index i of src 
[xli yli button] = ginput(l); 
20 if (button = 3) break; end 

d = sqrt((TOPO.Iocs(:,l) - xli).'^2 + (TOP0.1ocs(:,2) - yli).^2); 
[d,i] = min(d); 

xl = TOP0.1ocs(i,l); yl = TOP0.1ocs(i,2); 

25 

% find coords and index j of dst 
[x2i y2i]-ginput(l); 

d = sqrt((TOP0.1ocs(:,l) - x2i).^2 + (TOP0.1ocs(:,2) - y2i).^2); 
[d,j] = min(d); 

30 x2 = TOP0.1ocsC,l); y2 = TOP0.1ocsG,2); 

hold on; 

Ih = line([xl x2],[yl y2],'color';red'); 
35 cap = input('Enter capacity (in Mbps) > '); 

fprintf(l /About to create symetric %d Mbps link from node %d to node %d\n',cap,ij); 
doit = input('Enter Y to confirm, N to reject, and B to change bandwidth (Y)> Vs'); 
if (isempty(doit)) doit = '¥'; end 



if(doit=='n'|doit=='N') 
delete(lh); 
45 return; 
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end 



if(doit== V|doit = 'B') 

buf = sprintf('Enter capacity from %d to %d (in Mbps) > ',i,j); 
cap_i__to J = input(buf); 

buf = sprintf('Enter capacity from %d to %d (in Mbps) > 'j,i); 
capJ_to_i = input(buf); 



else 

1 0 cap_i_to J = cap; 

cap J_to_i ~ cap; 

end 

%build the link records 
1 5 clear linkab linkba; 

linkab.src = i; 
linkab,dst = j; 
linkab,bw = cap_i_toJ; 
20 linkab,handle = Ih; 

linkba.src = j; 
linkba.dst i; 
linkba.bw = capjjoj; 
25 linkba^handle = Ih; 



% now draw the actual link on the map 
delete(lh); 

30 Ih = drawlink(TOPO, linkab); 

% now store the link info 

TOPO.links =^ [TOPO.links ; linkab ; linkba ]; 

TOPOJinkarray = [TOPO.linkarray ; [ i j cap_i_toJ] ; [ j i capj_to_i ]]; 



end % of while loop 
assignin('caller',inputname(l),TOPO); 



function [C, portmap] ~ capacity(TOPO) 
45 % [C, portmap] = capacity(TOPO) 
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% portmap maps indices of C to elts of nodes(TOPO) 
% [node dir] where 

% node is index of elt in nodes(TOPO) 

% dir is 1 if data enters here, -1 if data leaves here 

5 

numnodes = length(TOPO. links) 2; 

C = zeros(nuinnodes,numnodes); 

10 curnode = 0; 

portmap = []; 

for i = l:length(T0P0.1inks) 
link = T0P0.1inks(i); 
cimiode - cumode + 1 ; 
1 5 portmap(cumode, :) = [link.src - 1 ] ; 

cumode = cumode + 1 ; 
.portmap(cumode,:) = [link.dst 1]; 

C(cumode-l,curnode) = link.bw; 

20 end 

c_node = 1; 
c_dir = 2; 

25 for i = 1 :length(TOPO.nodes) 

ins = fmd(portmap(:,c_node) ^= i & portmap(:,c_dir) == 1); 
outs = find(portmap(:,c_node) = i & portmap(:,c_dir) = -1); 

forj = ins 

30 for k = outs 

CO,k) = inf; 

end 

end 

end 

35 function [a, b, c] debug(t) 

update(t); 
fieldnames(t) 

40 

a = t.nodes 
b = t.locs 
c = t. links 

function dispIay(TOPO) 
45 % DISPLAY a topo object 
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% a link is a unidirectional, so the value is probably twice what you want 

fprintf('[TOPO object: %d nodes %d linksJW,... 

length(TOPO.nodes),length(TOP0.1inks)); 
5 function draw(TOPO) 

% draw(topo) 

% 

% draw the topology figure in a new window 

1 0 TOPO.cur_fig = figure; 

axis(TOPO.axis); 

axis equal; 

axis manual; 

box on; 
15 hold on; 

for i - 1 :length(TOPO.nodes) 

nm = plot(TOPO.nodes{i}.loc(l),TOPO.nodes{i}.loc(2);ob'); 
TOPO.nodes{i}.mark_handle = nm; 
20 if(isfield(TOPO,nodes{i};nameloc')) 

TOPO.nodes{i},nameloc(3) = text(TOPO.nodes{i}.nameloc(l),... 

TOPO.nodes{i} ,nameloc(2),TOPO.nodes{i} .name); 

end 

end 
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% yes, this draws the same link twice, fix it if it matters -dam 1 1/21 



T0P0.1inkarray-[]; 
for i = 1 :length(TOPO Jinks) 
30 T0P0.1inks(i).handle = drawlink(TOPO,TOPO,Hnks(i)); 

TOPO.linkarray = [TOPO.linkarray ; ... 

[ T0P0.1inks(i).src T0P0.1inks(i).dst TOPO.iinks(i).bw]]; 

end 

3 5 assignin('caller',inputname( 1 ),TOPO); 

function ex(t) 

t.nodes 

function labelnames(TOPO) 
40 % function labelnames(TOPO) 

% make it easy to label the nodes 

for i = l:length(TOPO.nodes) 

fprintfCPlace label for node %d "%s"\n',i,char(TOPO.nodes{i}.name)); 

45 
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origcolor = get(TOPO.nodes{i}.mark_handle/color'); 
set(TOPO.nodes{i}.mark_handle;color',[l 0 0]); 

if(isfield(TOPO.nodes{i};nameloc')) 

good_x = TOPO.nodes{i}.nameloc(l); 
good_y = TOPO.nodes{i}.nameloc(2); 

end 

th=[]; 
while (1) 

fprintf('Button 1 to (re)place text, Button 3 to accept\n'); 

[x,y,button] ginput(l); 

if (3 == button) break; end 

if (-isempty(th)) delete(th); end 

th = text(x,y,TOPO.nodes{i}.name); 

good_x = x; goodly = y; 

end 

TOPO. nodes {i}.nameloc = [good_x, goodly, th]; 
set(TOPO.nodes{i}.mark__handIe/color',origcolor); 

end 

assigninCcalIer',inputname(l),TOPO);function names(TOPO) 
% NAMES the list of names of the nodes in the topo 

fprintf('Node\t\tName\n'); 
for i = l:size(TOPO.names,l) 

fprintf('%d\t\t%s\n^i,TOPO.names{i}); 

end 

function [node] ^ nodes(TOPO) 

% function [node] = nodes(TOPO) 

% returns a cell array describing nodes in the TOPO 

node = TOPO.nodes; 

function [TOPO] - topo(TOPO) 

%[TOPO]=topo(TOPO) 

%% if input TOPO is 'init', create a new topology 

% 

% newtopo = topo('init'); 

% 

% else add new nodes to TOPO 
% 

% nodes is a array of structs, one per node 
% link is a array of structs, one per link 

% a link is a unidirectional item, so there are probably twice 
% as many links as you'd expect. 
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error('topo(TOPO) or topo("init") - not enough args'); 



if (ischar(TOPO) & TOPO 'init') 
clear TOPO 

TOPO.nodes = []; 
TOPO.links = []; 

TOPO.capacity = []; % now computed as needed 
TOPO.locs = [] ; % internal cache 
TOPO.linkarray = [] ; % internal cache 

f= figure; 
axis([0 75 0 50 ]); 
TOPO. axis = axis; 
TOPO.cur_fig = f; 
axis equal 
axis manual 
box on 
hold 

figure(TOPO.cur_fig); 



else 
end 



nodecount - length(TOPO.nodes); 

while (1) 

clear nodeinfo; 

fprintf(l,'\n\nHit Button 3 to stop\n\n'); 

[xybut] = ginput(l); 

if (but ==3) break; end 

X = floor(x); y = floor(y); 

nm = plot(x,y,'ob'); 

name = input('Enter name > '/s'); 

nodeinfo.loc = [ x y]; 
nodeinfo.markhandle = nm; 
nodeinfo. name = cellstr(name); 
nodecount = nodecount + 1 ; 
TOPO.nodes{nodecount} = nodeinfo; 



end 
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ifCtopo'-^= class(TOPO)) 

TOPO = class(TOPO;topo'); 

end 

5 

if (nargout =^ 0) 

assignm('cailer',inputname(l),TOPO); 

end 

function Ih = drawlink(TOPO, link) 
1 0 % assumes TOPO.linkarray is already valid, and draws the position of 

% link line based on the number of links already present in linkarray 

c_src = 1 
c_dst = 2 
15 c_bw==3 

i = link.src; 
j = link.dst; 

20 xl - TOPO,nodes{i}.loc(I) 

yl =TOPO.nodes{i}.loc(2): 
x2 = TOPO.nodes{j}.loc(l) 
y2 - TOPO,nodes{j}Joc(2) 

25 

if (isempty(TOPO-linkarray)) 
num_links = 0; 

else 

num_links = sum(T0P0.1inkarray(:,c_src) == i & T0P0.1inkarray(:,c_dst) j); 

30 end 

pattem=[0 1-12-2 3-3]*.3; 

if (abs(xl - x2) > abs(yl - y2)) 
35 delta_x = 0; 

delta_y = pattem(num Jinks + 1); 

else 

delta_x = pattem(num_links + 1); 
delta_y = 0; 

40 end 

Ih = line([xl x2] + delta_x, [yl y2] + delta_y, 'color', 'black'); 
function update(TOPO) 

45 



45 

clear TOPO.locs; 
fori=l:length(TOPO.nodes) 

TOPOJocs(i,:) = TOPO.nodes{i}.loc 

end 

5 

clear TOPO.linkarray; 

for i = l:length(TOPO Jinks) 

TOPO.linkarray = [TOPO.linkarray ; ... 

[ T0P0.1inks(i).src TOP0.1inks(i).dst T0P0.1inks(i).bw]]; 

10 end 

% these are here to be cut and pasted into other functions as needed 
% there doesn't seem to be a good way to pass them around in another fashion 
% (using assigning('caller\..) to force their definition sounds like asking 
1 5 % for trouble 'cause you'll overwrite another definition of them.,.) 

c_src = 1 ; 
c_dst = 2; 
Q c_J)w = 3; 

U 20 assigninCcaller'4nputname(l),T0P0); 



